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Abstract 

We demonstrate that recent advances in nanoscale thermal transport and temperature manipulation can 
be brought to bear on the problem of tailoring thermal radiation from compact emitters. We show that 
wavelength-scale composite bodies involving complicated arrangements of phase-change chalcogenide 
(GST) glasses and metals or semiconductors can exhibit large emissivities and partial directivities at mid- 
infrared wavelengths, a consequence of temperature localization within the GST. We consider multiple 
object topologies, including spherical, cylindrical, and mushroom-like composites, and show that partial di¬ 
rectivity follows from a complicated interplay between particle shape, material dispersion, and temperature 
localization. Our calculations exploit a recently developed fluctuating-volume current formulation of elec¬ 
tromagnetic fluctuations that rigorously captures radiation phenomena in structures with both temperature 
and dielectric inhomogeneities. 
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The ability to control thermal radiation over selective frequencies and angles through complex 
materials and nanostructured surfaces yj] has enabled umrecedented advances in important tech¬ 
nological area^ including remote temperature sensing IPQ, incoherent sources and energy¬ 

harvesting P-ItQ. Recent progress in the areas of temperature management and thermal transport 
at sub-micron scales can play a significant (and largely unexplored) role in the design of specially 
engineered radiative structures that combine both photonic and phononic design principles [8, 

In this letter, we describe a fluctuating-volume current (FVC) formulation of electromagnetic 
(EM) fluctuations that enables fast and accurate calculations of thermal radiation from complex 
structures with non-uniform temperature and dielectric properties. We demonstrate that when se¬ 
lectively heated, wavelength-scale composite bodies—complicated arrangements of phase-change 
materials and metals or semiconductors—can exhibit large temperature gradients and partially 
directed emission at infrared wavelengths. For instance, micron-scale chalcogenide (GST) hemi- 
spheroids coated with titanium or silicon-nitride shells [Fig. [U and resting on a low-index, trans¬ 
parent substrate can exhibit large emissivity and > 80% partial directivity—redirecting light away 
from the metallic or toward the semiconducting shell—when heated to 870 K by highly conduc¬ 
tive 2D materials at the GST-substrate interfaces. The interplay of geometry and temperature 
localization allows such composite infrared thermal antennas to not only enhance but also selec¬ 
tively emit and absorb light in specific directions. We show that other designer bodies, including 
mushroom-like particles and coated cylinders [Fig.|3l, can also exhibit large partial directivity, in 
contrast to situations involving homogeneous bodies or uniform temperature distributions which 
lead to nearly isotropic emission. Our predictions are based on accurate numerical solutions of the 
conductive heat equation and Maxwell’s equations, which not only incorporate material disper¬ 
sion but also account for the existence of thermal and dielectric gradients at the scale of the EM 
wavelength, where ray optical descriptions are inapplicable. 

Attempts to obtain unusual thermal radiation patterns have primarily relied on Bragg scat- 


terin g a nd related interference effects in nanostructured surfaces 


ings il 
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materials la 


metasurfaces Ill2l - ll9f] . multilayer structures mm, and sub-wavelength meta- 


faces 11 

iBS 


ij], including photonic grat- 


2414271]. Related ideas can also be found in the context of fluorescence emission. 


where directivity is often achieved by employing metallic objects (e.g . plasmonic antennas) to re¬ 


direct emission from individual dipolar emitters via gratings 11281 


1% or by 

31 


localizing fluorescent 

n, 

1150 • Matters become 


molecule(s) to within some region in the vicinity of an external scatterer y. 
complicated when the emission is coming from random sources distributed within a wavelength- 
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scale object, as is the case for thermal radiation, because the relative contribution of current sources 
to radiation in a particular direction is determined by both the shape and temperature distribution 
of the object. Optical antennas have recently been proposed as platforms for control and design 


of narrowband emitters Obi. 


3711 . though predictions of large directivity continue to be restricted 


to periodic structures. While there is increased focus on the study of light scattering from sub¬ 


wavelength partic 
color routing |l40|, 


es and microwave antennas (useful for radar detection 113 8ll . sensing 09tl . and 


4111 1. similar ideas have yet to be translated to the problem of thermal radiation 


a. 


Here, 


from compact, wavelength-scale objects, whose radiation is typically quasi-isotropic 
we show that temperature manipulation in composite particles could play an important role in the 
design of coherent thermal emitters. 

Temperature_gradients can arise near the interface of materials with highly disparate thermal 
]. Although often negligible at macroscopic scales 1420, recent experiments re¬ 


conductivities 


veal that the presence of thermal boundary resistance 1143 


4411 (including intrinsic and contact 


resistance 11450 ') in nanostructures together with large dissipation can enable temperature local¬ 
ization over small distances 11420 . Such temperature control has been recently investigated in the 
context of metallic nanospheres immersed in fluids [46], graphene transistors Hi , nanowire re¬ 


sistive heaters 11480 . AFM tips 11490 . and magnetic contacts 


high-symmetry structures, e.g. spheres 10] and planar films 123| 


. With the exception of a few 
1, however, calculations of ther¬ 


mal radiation from wavelength-scale bodies have been restricted to uniform-temperature operating 


conditions, exploiting Kirchoff’s law |123 


520 to obtain radiative emission via simple scattering 


calculations. The presence of temperature and dielectric inhomogeneities in objects with features 
at the scale of the thermal and EM wavelengths call for alternative theoretical descriptions. 

Formulation .— In what follows, we present a brief derivation of our FVC formulation of ther¬ 
mal radiation, with validations and details of its numerical implementation described in a separate 
manuscript 115.30 . Our starting point is the VIE formulation of EM scattering ll54ll . describing scat¬ 
tering of an incident, 6-component electric (E) and magnetic (M) field = (E; H) from a body 
described by a spatially varying 6x6 susceptibility tensor x(x). (Eor convenience, we omit the 
frequency cc dependence of material properties, currents, fields, and operators, and also define % 
to be the susceptibility relative to the background medium.) Given a 6-component electric (J) 
and magnetic (M) dipole source cr = (J; M), the incident field is obtained via a convolution (x) 
with the 6x6 homogeneous Green’s function (GF) of the ambient medium Efx. y), such that 
0inc = r X 0 - = J d^yr(x, y)a(y). Exploiting the volume equivalence principle ll54ll . the un- 
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Figure 1. Schematic of representative composite bodies comprising GST (blue) hemispheroids coated with 
Ti (green), AZO (red), or Si 3 N 4 (orange) shells, and resting on a low-index, transparent substrate in contact 
with a heat reservoir at 300 K. The GST is heated from below by a conductive 2D material (e.g. a carbon- 
nanotube wall or graphene sheet), leading to temperature gradients within the structure. The presence of 
boundary resistance at material interfaces is captured by effective (intrinsic and contact) thermal resistances 

^th- 

known scattered fields (/)sca = T ^ can also be expressed via convolutions with F, except that 
here ^ ^^e the (unknown) bound currents in the body, related to the total field inside the 

body (/) = (/)inc + </>sca through x- Writing Maxwell’s equations in terms of the incident and induced 
currents, 


i *0 = -*^x(r * cr), 


( 1 ) 


one obtains ^ in terms of the incident source a. 

Consider a Galerkin discretization of Eq. [U via expansions of the current sources cr(x) = 
En Sn&n(x) and ^(x) = En Xnbni'x.) in a convenient, orthonormal basis {bn} of N 6-component 
vectors, with vector coefficients s and x, respectively. The resulting matrix expression has the 
form X -\- s = Ws, where = (^m, bn + «<^x(r 'k bn)) is known as the VIE matrix and 

(,) denotes the standard conjugated inner product. Poynting’s theorem implies that the far-field 
radiation flux $ = ^ Re J d^x (E* x H) due to a. 


Rer0 = - J Re(^ + a)*T x (^ + a) 


= —+ -s)* sym G{x -|- s) = —-s*W* sym GWs 

Zd Zi 
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where s*s and G are both N x N matrices and Gmn = (^m, F * ) are the elements of the 

so-called “Green” matrix. Thermal radiation from such a body follows from the cumulative flux 
contributions of a collection of incoherent sources distributed throughout its volume 115311 . obtained 
by a thermodynamic, ensemble-average $ = ($ 0 ^) over all a and polarizations. It follows that the 
total radiation is given by: 

$ = -ixr [GW* sym GW], (2) 

where we have defined the current-current correlation matrix C, whose elements Cmn = {sm^n*) = 

f f d^xd^v bZj(x)(a(x)a*( y)) br,(v). The correlation functions satisfy a well-known fluctuation- 
dissipation theorem (FDT) 1155n . (ai(x,u)a*(y,Lo)) = Im x(x, cj)0(x, cj)S(x — y)Sij, relating 
current fluctuations to the dissipative ^ Imx and thermodynamic properties of the underlying 
materials. Here, 0(x, a;) = — 1) is the Planck distribution at the local temperature 

T(x) 115fin . In addition to the total flux, it is also desirable to obtain the angular radiation pattern 
from bodies, which can be straightforwardly obtained by introducing the far-field Green matrix 
= {bjn, QF^ -k bn), based on the 3x6 GF F^(x, y) which maps currents to far-field electric 
fields and the 3 x 3 tensor Q mapping Cartesian to spherical coordinates (azimuthal and polar 
angles) [? ]. Following a similar derivation, the angular radiation flux in a given direction will be 
given by: 


U = 


k^Z 

2(47r)^ 


Tr [CW(G^^Gl)W ], 


(3) 


where k — uj/c and Z — \/iiq/sq is the wave impedance of the background medium. Equation [3] 
can be employed to calculate emission from arbitrarily shaped bodies with spatially varying di¬ 


electric and temperature properties—unli 


ce previous scattering-matrix and surface-integral equa¬ 


tion formulations of thermal radiation ll56h . the FVC scattering unknowns are volumetric currents. 


Furthermore, when considered as a numerical method as we do below, the corresponding basis 
functions can be chosen to be localized “mesh” elements 1^ . allowing resolution to be employed 
where needed l53n . 

Results .— In what follows, we explore radiation from composite bodies comprised of chalco- 
genide Ge 2 Sb 2 Te 5 (GST) alloys and metals or semiconductors. To begin with, we consider micron- 
scale GST hemispheroids coated with either titanium (Ti) or silicon-nitride (Si 3 N 4 ) shells, de¬ 
picted in Fig. [U and described in detail in [SM]. The structures rest on a low-index (e 1) 
transparent substrate which not only provides mechanical support but also a means of dissipat¬ 
ing heat away from the structure; the bottom of the substrate is assumed to be in contact with 
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Figure 2. Spectral emissivity e{Ljj) = $(6c;)/$bb(<^) (blue dots) from composite bodies comprising 
Ge 2 Sb 2 Te 5 (GST) hemispheroids coated with either (a) Si 3 N 4 or (b) Ti shells, under heating scenarios 
(ii) and (hi), respectively (see text). The dashed blue line in (a) shows e for the Si 3 N 4 structure under 
heating scenario (iii). Both structures rest on a low-index substrate (shown schematically in Fig. [T]) and 
are heated from the GST-substrate interface by a 2D thin-film conductor, reaching an interface temperature 
Tgst = 870 K; the temperature profiles r(x) are shown as insets, e is defined as the ratio of the thermal 
flux $(6t;) from each body normalized to the flux $bb = 4 ^(^/c)^©(<^, T) from a corresponding black 
body of same area A and uniform temperature T — 870 K (green lines, arbitrary units). Also shown are the 
partial directivities r]± = (red line), defined as the ratios of the outgoing flux into the upper/lower 

hemisphere = 27r 9) to the total flux, where 9 is defined in Fig. [T] (c) Angular radia¬ 

tion intensity U{9) for the structures in (a) (solid red line) or (b) (solid blue line) as well as under heating 
scenario (v), corresponding to uniform temperature distributions (dashed lines), (d) Total (frequency in¬ 
tegrated) downward partial directivity P_ = H^/H as a function of Tqst, where i7_ = dcj^-(cj) 
and H — duo $(cc;), for the Ti structure under different heating conditions, corresponding to multiple 
degrees of temperature localization in the GST (see text). 


a 300 K heat reservoir while surfaces exposed to vacuum satisfy adiabatic boundary conditions 
(VT • fi = 0). When heated by a highly conductive 2D material (e.g. carbon nanotube wall 
or graphene sheet) at the GST-substrate interface, such a structure can exhibit large temperature 
gradients within the core, a consequence of boundary resistance between the various interfaces 


and rapid heat dissipation in the highly conductive shells 0431 bSL l59l]. To model the correspond¬ 


ing steady-state temperature distribution T(x), we solve the heat-conduction equation via COA - 
SOL |3], including the full temperature-dependent thermal conductivity k{T) of the GST 
Note that even at large temperatures, «;gst(800 K) <C KTi,Si3N4(300 K) > 20 IL/mK. The ex¬ 
istence of (intrinsic and contact) boundary resistance at this scale is taken into account by the 
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Figure 3. Angular radiation intensity U{9) for a variety of heterogeneous bodies—composite GST(blue), 
AZO (red), and Si 3 N 4 (orange) ellipsoids, spheres, mushrooms and cylinders—at selected frequencies u. 
The temperature of the GST is held at 870 K while that of other materials is fixed at 300 K. e and ri± denote 
the emissivity and partial directivities in the upper/lower hemispheres, defined (along with 6) in Fig. [T] and 
in the text. Polar (3D) plots are normalized by the total flux. 


introduction of effective resistances Rsh\c, ^h|su, and Rsh\su, at the interfaces between shell-GST, 
heater-substrate, and shell-substrate, respeetively [SM]. Figure |2td) shows T’(x) throughout the 
Ti structure when the GST-substrate interfaee is heated to Tqst = 870 K (approaching the GST 
melting temperature ll62n ). and under various operating conditions. Specifically, we consider 
^sh|su = 10“®m^FF/K and i?sh|c = -Rh|su = Rth, where {(i), (ii), (iii)} correspond to typical 
values of Rth = {0.5,1,2} x 10“^m^FF/K while (iv) Rth = oo and (v) Rth = 0 deseribe extreme 
situations involving either perfect temperature localization in the GST or uniform temperature 
throughout the strueture, respectively. 


Given T’(x) and the dieleetric properties of the bodies 1163 - 


6611 ■ the flux can be obtained 


via Eq. [3l Due to temperature gradients and phase transitions in the GST 116 lh . its dielectric re¬ 
sponse £gst(7’(x), C(j) eonsists of continuously varying rather than piece-wise eonstant regions 
[SM]; our FVC method, however, can handle arbitrarily varying £(x) and T’(x). The ehoiee of 
materials, shapes, and dimensions of the hemispheroids ensure the existenee of geometric reso¬ 
nances near the thermal wavelength At ~ 5.8/um, corresponding to the peak of the blackbody 
spectrum at 870 K. In this wavelength regime, £gst ~ 30 -I- ll63l . 16411 . £Ti ~ —100 -I- 80) llOfill . 
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and £si 3 N 4 ~ 5 + O.li ll66h . enabling significant Purcell enhancement and emission ll67h . Fig¬ 


ure |2] shows the emissivity and partial directivity of the (a) Si 3 N 4 and (b) Ti structures, along 
with their corresponding T(x) (insets) under heating scenarios (ii) or (iii), respectively, and for 
Tgst = 870 K. The emissivity e{uj) = $(a ;)/$bb (<^) (blue dots) is defined as the ratio of the ther¬ 
mal flux $(c(j) from each object to that of a blackbody $bb(<^) = T) of the same 

surface area A and uniform T = 870 K (green lines); the partial directivity r]± = $± 7 $ (red lines) 
is defined as the ratio of the flux into the upper/lower hemisphere 

to the total flux $, where 6 is defined with respect to the +z axis [Fig. [H. Note that although 
e exhibits multiple peaks, its magnitude (e < 0.2) is limited by material losses (Ime < Ree) 
in this frequency range 11671] : larger e can likely be optained with further design and/or material 


combinations. 


We find that ij increases sharply as the system transitions from quasistatic to wavelength-scale 
behavior (in contrast to e which exhibits gradual variations, except near a resonance). At small 
LoLjc <C 1, the emission is highly quasi-isotropic (as expected from a randomly polarized dipolar 
emitter lll2l]i. becoming increasingly asymmetric as uL/c > 1. Essentially, with the help of the 


curvature [l68ll . the Si 3 N 4 and Ti shells redirect radiation upwards or downwards, enabling strong 


coherent interference between the radiated and scattered fields of dipole emitters within the GST, 
making the design of the temperature profile an essential ingredient for achieving large ry. Fig¬ 
ure l2tc) shows the angular radiation intensity U(6) of the Si 3 N 4 (red) and Ti (blue) structures 
at selected frequencies and under two of the above-mentioned heating conditions, corresponding 
to either (ii) partial temperature localization in the GST (solid lines) or (v) uniform temperature 
throughout the bodies (dashed lines). The dramatically different radiation patterns and signif¬ 
icantly smaller ry under (v) belie the fact that dipole emitters inside the GST contribute larger 
partial directivity compared to those in the shell, which tend to radiate quasi-isotropically and 
dominate e. 


To illustrate the non-negligible impact of T’(x) on the total radiation of the bodies. Fig. |2td) 
shows the frequency-integrated, downward partial directivity P_ = H_/H of the Ti structure un¬ 
der different Tqst and heating scenarios, where i7_ = du 'h_(a;) and H = duo $(a;). As 
expected, P_ grows with increasing temperature localization in the GST, and remains almost con¬ 
stant fa 0.55 under uniform temperature conditions. Such an increase in partial directivity, how¬ 
ever, comes at the expense of decreasing e (not shown) due to the increasingly dominant role of 
larger frequencies. At large uj or for large bodies (where ray optics becomes valid), material losses 
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severely diminish e. Not surprisingly, the design criteria of such wavelength-scale emitters differs 
significantly from that of large-scale bodies (where Kirchoff’s law is valid 11691] 1. For instance, 
while larger rj can be obtained in the ray-optics limit by increasing the shell thickness of each 
structure relative to the GST dimensions (thereby enhancing extraction/reflections of radiation 
from the core), optimal r] at a fixed frequency ujL/c ~ 1 occur at specific shell thicknesses, deter¬ 
mined by the shape and materials of the bodies. Planar structures can also yield highly directional 
emission when subject to inhomogeneous temperature distributions ll23n . but require significantly 
larger boundary resistance (heating power) and offer limited degrees of freedom for controlling 
emission. Compared to large-scale or planar radiators, wavelength-scale composite bodies not 
only provide a high degree of temperature tunability, but also enable simultaneous enhancement 
in rj and e, even potentially exceeding the ray-optical, blackbody limit ifiTH . 

Figure [3] shows the radiation pattern of other heterogeneous structures (ellipsoids, spheres, 
mushroom-like particles, and cylindrical composites), depicted schematically in the figure with 
blue/red/orange denoting GST/AZO/Si 3 N 4 materials. Their shapes and dimensions are detailed 
in [SM]. For simplicity, we consider emission at selected frequencies and under heating scenario 
(iv), corresponding to perfect temperature localization in the GST. As expected, the design criteria 
for achieving large rj differs depending on the choice of materials, with GST-Si 3 N 4 composites 
favoring large-curvature prolate bodies llfisll and GST-AZO composites favoring oblate structures 
that provide higher reflections. 

Concluding remarks .— The predictions above provide proof of principle that combining con¬ 
ductive and radiative design principles in wavelength-scale structures can lead to unusual thermal 
radiative effects. Together with our FVC formalism, they motivate the need for rigorous theoretical 
calculations of thermal emission that account for existence of temperature and dielectric gradients 
in micron-scale, structured surfaces, an issue that is especially relevant to thermal metrology H. 
The FVC framework not only enables fast and accurate calculations, but also for techniques from 
microwave antenna design and related fields to be carried over over to problems involving in¬ 
frared thermal radiation. Although the focus of this work is on thermal radiation, similar ideas 
and techniques are applicable to problems involving fluorescence or spontaneous emission where, 
rather than controlling the temperature profile, it is possible to localize and control the sources of 




emission via doping 117 If] or judicious choice of incident laser light 

We are grateful to Bhavin Shastri for very helpful comments. This work was supported in part 
by the National Science Foundation under Grant No. DMR-1454836. 
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Supplemental Materials: Temperature control of thermal radiation from 

heterogeneous bodies 

Below, we provide details of the geometric and material properties of the composite bodies 
described in the main text, along with addition discussion of the parameters and assumptions of 
the heating schemes leading to temperature gradients. 


1. GEOMETRY AND MATERIAL PARAMETERS 



Figure S1. Shape parameters describing the Si 3 N 4 (left) and Ti (right) composite hemispheroids explored 
in Fig. 2 of the main text. Blue, orange, and green colors denote GST, Si 3 N 4 , Ti, respectively. 


Figure ISH provides schematics of the Si 3 N 4 (left) and Ti (right) hemispheroid composites ex¬ 
plored in Fig. 2 of the text, along with the corresponding geometrical parameters. In particular, 
the Si 3 N 4 and Ti shells have long (short) semi-axes 2.6(1.9)/rm and 1.2(0.46)/rm, respectively; 
the chalcogenide (GST) cores have long (short) semi-axes of 0.88(0.63)/um and 0.83(0.33)/um, 
respectively. As noted in the text, these values are chosen in order to enhance the emissivity and 
partial directivity of the structures. Figure [S4] provides the size and dimensions of the various 
geometries explored in Fig. 3 of the text. 

The Ge 2 Sb 2 Te 5 alloy is a phase-change chalcogenide glass that exhibits a large thermo-optic 
effect 1^ and three possible (amorphous, cubic, and hexagonal) phases corresponding to transi¬ 
tion temperatures of 438 K (separating the amorphous and cubic phases) and 623 K (separating 
the cubic and hexagonal phases) 11581 . l64l . Because there are yet no experimental characteriza¬ 
tions or semi-analytical models of the dielectric dispersion £GST(<^,r) of the GST from 300 K 
to its melting point 870 K |^, we instead model the dispersion via a simple linear-interpolated 
fit of available experimental data at five different temperatures (spanning amorphous, cubic and 
hexagonal phases) ll63L 1641]. Figure [S2l shows both the real (red solid line) and imaginary (blue 
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Figure S2. Real (solid red) and imaginary (dashed blue) permittivity ecsT {T, A) of a bulk GST glass at 
A = 5.8/Lim, obtained via simple linear interpolation of experimental data at multiple temperatures (circles). 


dashed line) parts of sgst at a single wavelength A = 5.8;um over this temperature range (with 
circles denoting experimental data). Together with the temperature profiles of the structures T(x) 
and dispersion relations of Ti B, Si 3 N 4 B, and AZO B, this provides all of the informa¬ 
tion needed to perform the calculations of thermal radiation from the bodies of Fig. 2 in the main 
text. On the other hand, Fig. 3 of the main text explores single-frequency radiation from bodies 
with piece-wise constant temperature profiles (constant T = 870 K in the GST and 300 K in the 
remaining regions), which allows us to employ typical permittivity values for these materials at 
mid-infrared wavelengths; specifically, we choose £gst = 30 -|- lOi [^], s^i = —100 + 80i 1^ 1. 
£si 3 N 4 = 5 -h 0.1) IB? and £azo = “25 -|-15), corresponding to a doping density ~ lwt% 1741. 


II. TEMPERATURE GRADIENTS 


Interfaces play a crucial role in nanoscale thermal transport. For instance, they enable thermal 
boundary resistance (TBR) to radically alter the surrounding temperature distribution ll43l-l46ll. 
leading to small-scale thermal discontinuities across the interface. TBR consists of both con¬ 
tact and intrinsic “Kapitza” resistance, with the former arising from poor mechanical connection 
between materials (due to surface roughness) and the latter from acoustic mismatch between ma¬ 
terials (and hence persisting even under perfect-contact situations) ll45n . Typical values of intrinsic 


resistance at room temperature are on the order of 10 ® ~ 10 'm^lU/K ll45h . whereas those 
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Figure S3. Temperature distribution r(x) of the Ti (left) and Si 3 N 4 (right) hemispheroid composites de¬ 
scribed in Fig. 1. Both structures rest on a Si02 substrate (thickness 0.3nm and radius = 1.5x shell radius) 
whose bottom surface is in contact with a heat reservoir at 300 K. All other surfaces are exposed to vacuum 
and therefore satisfy adiabatic boundary conditions; material interfaces on the other hand are subject to 
thermal boundary resistance in accordance with operating condition (iii) and (ii) described in the text, for 
the left and right body, respectively. 


arising from contact resistance vary depending on the surface and thermophysical properties of 
the intervening medium. In our setup (described schematically in Fig. 1 of the main text), there 
four interfaces at which TBR can arise. These are denoted and described by the resistances Rsh\c, 
^h|su> ^sh|su> and i?h|c 5 of the shell-GST, heater-substrate, shell-substrate, and heater-GST in¬ 
terfaces, respectively. Note that the thermal resistance associated with graphene can be made 
extremely small IITSH and hence in our calculations, we assume negligible i?h|c = 0. In order to 
obtain large temperature gradients, it is important to operate with materials that can dissipate heat 


away from the shells rapidly hence, we assume small shell-substrate interface resistances 
Rsh\su = 10“®m^FF/K. For simplicity, we consider conditions under which the interface resis¬ 
tances Rsh\c = f^h|su = Rth are equal and obtain various degrees of temperature localization by 
varying Rth, with Rth = oo leading to perfect temperature-localization and Rth = 0 leading to 
uniform temperature distributions. In particular, we consider five different operating conditions, 
corresponding to realistic values of (i) Rth = 0.5 x 10“^m^FF/K, (ii) Rth = 10“^m^FF/K, and 
(iii) Rth — 2 X 10“^m^FF/K and extreme (unrealistic) values of (iv) Rth = oo, and (v) Rth — 0. 


In order to solve the heat-conduction equation to obtain the steady-state temperature distri¬ 
bution T’(x), one must also specify the boundary conditions associated with vacuum-material 
interfaces, which we assume to be adiabatic (VT • ri = 0), corresponding to negligible conduc¬ 
tion, convection, and radiative-heat dissipation through air. The substrate is chosen to be a 0.3 /im 
thick Si02 film in contact with a heat reservoir at 300 K through the bottom interface, leading to 
large heat dissipation away from the shell and hence large temperature localization in the GST 
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Table I. The thermal properties used in our COMSOL simulation 



k(W/m/K) 

C(J/kg/K) 

Ti 

21.9 

523 

Si3N4 
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Figure S4. Parameter descriptions for the bodies associated with Fig. 3 of the main text. 


with decreasing substrate thickness. We choose the substrate lateral (cylindrical) dimensions to be 
large enough to remove large thermal diffusion away from the GST. Figure [S3] shows the temper¬ 
ature distribution of both Ti (left) and Si 3 N 4 (right) structures under the operating condition (iii) 
and (ii), respectively, assuming the material conductivities and heat capacities given in Table I. As 
shown, the temperature in the substrate and shells is almost uniform and close to 300 K, thanks 
to the presence of boundary resistance between the heater and substrate (which bars heat from 
flowing into the substrate) along with the high thermal conductivities of Ti and Si 3 N 4 , which act 
to dissipate heat away from the GST. 
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